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Abstract 

The neutrino-nucleus cross section and the muon capture rate are discussed within a simple 
formalism which facilitates the nuclear structure calculations. The corresponding formulae only 
depend on four types of nuclear matrix elements, which are currently used in the nuclear beta 
decay. We have also considered the non-locality effects arising from the velocity-dependent terms 
in the hadronic current. We show that for both observables in ^■^C the higher order relativistic 
corrections are of the order of ~ 5% only, and therefore do not play a significant role. As nuclear 
model framework we use the projected QRPA (PQRPA) and show that the number projection plays 
a crucial role in removing the degeneracy between the proton-neutron two quasiparticle states at 
the level of the mean field. Comparison is done with both the experimental data and the previous 
shell model calculations. Possible consequences of the present study on the determination of the 
i^fi neutrino oscillation probability are briefiy addressed. 



PACS number: 23.40.-s, 23.40.Bw, 23.40.Hc, 25.30.Pt, 21.60.Jz 
Keywords: neutrino-nucleus cross section, muon capture, projected QRPA 



1 



I. INTRODUCTION 



The semileptonic weak interactions with nuclei include a rich variety of processes, such as 
the neutrino (antineutrino) scattering, charged lepton capture, decays, etc, and we have at 
our disposal the results of more than a half-century of beautiful experimental and theoretical 
work. At present their study is mainly aimed to inquire about possible exotic properties 
of the neutrino associated with its oscillations and massiveness by means of exclusive and 
inclusive scattering processes on nuclei, which are often used as neutrino detectors. An 
example is given by the recent experiments performed by both the LSND and the 

KARMEN j3| Collaborations, looking for — * z/g and oscillations with neutrinos 

produced by accelerators. When the Vf. come from decay-at-rest (DAR) of /i"*" the flux 
contains neutrinos with a maximum energy energy of 50 MeV and can be detected through 
both exclusive and inclusive Vf. + ^'^C ^'^N + e~ reactions |i^J. In the case of coming from 
the decay-in-flight (DIF) of vr"*", the neutrino flux extends over the range (0,300) MeV, and 
the v^j, —>■ Ue appearance mode is looked for ejgjerimentally through the reaction u^ + ^'^C — *• 
^^A^ + which also has been measured [5]. On the other hand we don't have at our 
disposal experimental information on the z/g-reaction in the DIF energy range, which is 
necessary in the evaluation of the oscillation probabilities. Therefore, it is imperative to 
develop nuclear models, capable of reproducing the Ue-DAR and z/^-DIF data, to be used to 
predict reliable values for the + ^"^C ^"^N + e~ cross section in the DIF energy range, 
and to calibrate forthwith the z/^ Vf. appearance probability. Needless to say, in addition 
and for consistence, the implemented model should also describe properly the well-known 
^^A^ ^"^C + z/g -|- e+ /^"'"-decay and the /i~ + ^^C ^ z/^ + ^^B muon-capture modes. 

The weak interaction formalism most frequently employed in the literature is that of 
Donnelly and Walecka where the nuclear form factors are classifled as Coulomb {M), 

longitudinal (£), transverse electric (T^') and transverse magnetic (T™'*^), in close analogy 
with the electromagnetic transitions. They in turn depend on seven nuclear matrix elements, 
denoted as: Mf,Af,A' 

J ) Sj^i S'j ) ^" J ^-iid QY . However, in studying neutrino induced 
reactions 0, 0] it is sometimes preferred to employ the formulation done by Kuramoto et 
a/.^^, mainly because of its simplicity. The later formalism does not include the velocity 
dependent terms in the hadronic current, and therefore the reaction cross section only de- 
pends on three nuclear form factors, denoted as |(/|cr|i)p and A. What is more. 
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the formalism of Kuramoto et a/.Q] does not include the muon capture rates. Therefore, 
to describe simultaneously the neutrino-nucleus reactions and /i-capture processes it is nec- 
essary to resort to additional theoretical developments, such as those of Luyten et a/. 



and Auerbach and Klein where one uses the matrix elements My, M\ and Afp that are 
related to the former ones in a non-trivial way. 

Quite recently, we have carried out a multipole expansion of the V — A hadronic current, 
similar to the one .ed Ba..„o . S do.He beta deca. 

expressing all above mentioned observables in terms of the vector (V^) and axial vector 
{A) nuclear form factors M.v and M.ai with M = — 1,0, +1 Such a classification is 

closely related to the V — A structure of the weak current and to the currently used nuclear 
/9-decay formalisms, where the nuclear moments are expressed in a way that the forbiddeness 
of the transitions is easily recognized. This in no way means that we have developed a new 
theoretical framework; the final results can be found in one form or another in the literature. 
The main difference stems from the fact that we use the Racah algebra from the start, 
employing the spherical spatial coordinates instead of the cartesian ones. Concomitantly, 
we also express the lepton trace in spherical coordinates, as done for instance in the book of 
Supek . Here we present a few more details than in Ref . 3| on the procedure which we 
have followed and also consider the first order nonlocal corrections, which give rise to the 
additional matrix elements Aiyi and Ma'- 

In Ref. we have analyzed the inconveniences that appear in applying RPA-like models 
to describe the nuclear structure of the {^^B, ^^C, ^^A^} triad. We have established that the 
projected quasiparticle RPA (PQRPA) was the proper approach to treat both the short 
range pairing and the long range RPA correlations. More details are given in the present 
work, putting special emphasis on the differences between the projected and the usual QRPA 
approximations. We also compare our PQRPA cross sections with recent shell model (SM) 
calculations Q, [l^; and analyze the reliability of the calculated energy dependence of the 
Ue + ^"^N + e~ cross-section in the DIF energy region to be used in the evaluation of 

the oscillation probabilities. 

The work will be organized as follows. In Section 11 we explain our multipole decomposi- 
tion of the weak current and present the corresponding formula for the neutrino-nucleus cross 
section and the /i-capture. In Section III we briefly overview the projected BCS (PBCS) 
and the PQRPA equations. The main reason for the last is to emphasize the relationship 
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with the usual QRPA approximation, which is difficult to find in the existing literature. 
Finally, in Section IV we show the numerical results which come from the PQRPA model 
and compare them with the most recent shell model calculations. 



II. FORMALISM FOR THE WEAK INTERACTING PROCESSES 

The weak Hamiltonian is expressed in the form 0, Q, ^ , 18 1 

H^{v) = ^JJ^e-"-\ (2.1) 
where G = (3.04545 ± 0.00006) xlO~^^ is the Fermi coupling constant (in natural units). 



2M m£ 



(J,zj0) (2.2) 



is the hadronic current operator {4^, and 

/a = -iUs^{Yi,El,)'-f^{l + '-f^)usS^,Ey) = {\il^li) (2.3) 

is the plane wave approximation for the matrix element of the leptonic current in the case 
of neutrino reactions. Here = 1,2,3,4, and Walecka's notation [3] with the Euclidean 
metric for quadrivectors is employed, i.e. x = {x, X4 = ixfn}. The only difference is that we 
substitute his indices (0,3) by our indices (0,0), which means that we use the index for 
the temporal component and the index for the third spherical component. The quantity 

k = Pi-Pf = {k,ikfi,} (2.4) 

is the momentum transfer {Pi and Pf are momenta of the initial and final nucleon (nu- 
cleus)), M is the nucleon mass, me is the mass of the charged lepton, and gv, Qa, Qm and 
Qp are, respectively, the vector, axial-vector, weak-magnetism and pseudoscalar effective 
dimensionless coupling constants. Their numerical values are 

Qv = I; gA = 1-26; 

2Mm, 

Qm = i^p- nn = 3.70; gp = S'at^- r- (2.5) 

+ m^ 

These estimates for ^f^f and Qp come from the conserved vector current (CVC) hypothesis, 
and from the partially conserved axial vector current (PCAC) hypothesis, res pec tively. In 



the numerical calculation we will use an effective axial-vector coupling 



1 QQ, 



2l| 



The finite nuclear size (FNS) effect is incorporated via tlie dipole form factor with a cutoff 
A = 850 MeV, i.e. as: 



A2 



9^9 



A2 + A;2 ; 

To use p.lj) with the non-relativistic nuclear wave functions, the Foldy-Wouthuysen trans- 
formation has to be performed on the hadronic current (|2.2|) . When the velocity dependent 
terms are included this yields: 

Jut = 9v + {9a + 9pi)(^ ■ k - 9aO- ■ V, 
J = -g^o- - ig^a- x k - ^^k + gp^icr ■ k)k + g^v, 

(2.6) 

where v = — iV/M is the velocity operator, acting on the nuclear wave functions. The 
following short notation has been introduced: 



(2.7) 



2Mm/ ^"^^ ^ 2Mm/ 
where k = |k|. 

In performing the multipole expansion of the nuclear operator 

0^ = {0,iO^) = J^e-'''-\ (2.8) 
the momentum k is taken to be along the z axis, i.e. 



^rV47r(2L + l)ji(/tr)llo(f), (2.9) 

L 



and one gets 



00 = Y.^~y^7^{2J+l)jj{Kr)Yjo{r)J9, 
J 



Om = ^iiMV4vr(2J + l)jL(/s:r)[r,.(f)® J]^^. (2.10) 

JL 



The geometrical factors 

/ 



L 1 J 
-M M 



(2.11) 
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TABLE I: Values of the geometrical factors Fjlm- 



M 



-1 



J + 1 



J 



J- 1 



J + 1 



J 



J- 1 



J + 1 



FjLM 



V 



J+1 

2 J+1 



\/ 2J+1 
a/2 V 2J+1 



1 



J+1 

A^V 2J+1 



1 / J 



V2 V 2J+1 



-1 



J 



1 

a/2 



J- 1 



J_ ./ J+1 
V2 V 2J+1 



which will be seldom used in our work, fulfill the sum rule 

FjlmFjlM' — ^MM', 

L 

and their explicit values are listed in Table I. 
After some Racah algebra we find 



(2.12) 



J 



with 



(2.13) 



O0( J) = i [^vYjo(«r) - c/AYjo(«r, <r • v)] + Vivf (^^ + ^pj ^ i ^FjloSjloM 

(2.14) 



L=J±1 



and 



Om(J) = ViTT^i ^-FjLM [-5'ASjLM(«r) + g'vPjLM(«r) - i(-)-^g'v-^jLoYjM(«r) 



(2.15) 



TABLE II: Elemental operators and their parities. 



Operator 


Parity 


YjM(Kr) = jj{Kr)YjM{r) 


i-y 


SjLMinr) = jiinr) [Yiir) ® a]jj^j 




PjLMinr) = jLiKr)[YL{r) v]jm 




YjMinr, <T • v) = jj{Kr)YjM{r)icr ■ v) 





where we have introduced the operators listed in Table ITTl and the coefficients 

/ 



G 



JLI 



6(2L + 1)(2/+1^ 



1 1 1 
I L J 



I 1 L 




which obey the sum rule 

FjlmGjli = —MFjiM- 

L 

Using (j2.12j) and (j2.17|) we can rewrite the spherical components of Om{J) as 

-gj, + Mg^ + gp25Mo)SjLM{fir) + gyPjLMif^r) 



Om{J) = V4^Y.^~'^FjLM 
L 

- V4Tri~-^g^SMo'^JoM- 

For the neutrino-nucleus reaction the momentum transfer is k = p 
{pi,iEi} and q^, = {qy,iE^}, and the corresponding cross section reads 

a{Ee,Jf) = + f d{cose)%{K,Jf), 

2,71 J -I 



where F{Z + 1, E() is the Fermi function, 6* = q^, ■ p^, and 



(2.16) 



(2.17) 



(2.18) 
q^, with pi = 



(2.19) 



(2.20) 



with I JjMj) and \JfMf) being the nuclear initial and final state vectors. 
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The weak hamiltonian matrix element reads 



{JfMf\H^\JM) = -^(^Ja^ 

where 



(2.21) 



a = {JfMf\0^\JM), 
and la is the leptonic current defined in ()2.3p . Thus 

Jf) = ^aO^Ca^, 

where, the lepton trace £q,^, when expressed in cartesian spatial coordinates, reads 



(2.22) 



(2.23) 



SlSv 



1 



EgEy 



(2.24) 



the positive (negative) sign standing for neutrino (antineutrino) scattering. 



It is convenient to follow Ref. 



15[ and express the spatial parts of O and C in spherical 



coordinates (M, M' = 0, —1, 1). In this way one might write 



with 



2 7-1-1 ^ 



C©!^ + X! ^mO\pCmm' — 23? I O0 ^(— )*'^C_jv/£0jv/ 

MM' \ M y 



(2.25) 



Cm = -C, 



44 



p- q 



(2.26) 



4M 



EiE^ 



[EeqM + EypM =F »i(q x p' 



M 



(2.27) 



MM' 



— ^MM' + 



1 



(?mPm' + qM'P*M - hiM'P ■ q) 

q*M-M' P*M-M' 



ElEy 

1 1 1 

■M M' M - M' 



E, 



Ep 



(2.28) 



From the Wigner-Eckart theorem we also get: 



J f J Ji 



J \ -Mf M,, 



(2.29) 
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and 



J 

Using now the orthogonality condition 
/ 



{ Jf J J^\ 

\-Mf MMij 



{Jf\\OM{J)U). 



(2.30) 



E 

MiMf 



Jf J Ji \ I J f J' Ji 



f 



1 



y-Mf M Mi j \-Mf M' M, / (2^ + 1) 



SjJ'^MM' 



(2.31) 



one obtains 



Jf) 



2J^ + 



\{jf\mjm)?c„+Y.\{jf\\OM(jm)\^CMM 

M 

-2^{{Jf\mj)Unjf\\Oo{J)U)C9o)], (2.32) 



where, from (imi and (ITT^ . 



{jf\mj)U) 

(^/I|Om(J)||J.) 



2Ji + l 
2Ji + 1 



gyMyiJ) - g^MA'{J) + (^a + ^pJ MI{J) 
{-g^ + Mg^ + g^,6Mo)Mf{J) + g^M^^J) - gjMoMv{J) 



(2.33) 



An 



with the nuclear matrix elements defined as 
Mv{J) 

Ml',{J) 



I An 
2J~+ 



——{Jf\\Yj{Kr)U), 



jj:^~''FjLM{Jf\\SjL{^r)U), 



(2.34) 



Air 



^j—^{J^\\Yj(nr,a.v)\\J.^, 



An 



2Ji 



-Y.^-'^FjLM{Jf\\PjL{^v)U). 

-L T 



The explicit expressions for Ctj^M and Cmw that appear in ()2.32|) are jl^ : 

/tin T)o \ 

= 

-^00 = 



Eu El, 
, 2goPo - P ■ q 



-±i,±i 



E(,Ey 

Ev El 



EiEy 



(2.35) 
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with 



^0 

Po 



■ q 

k ■ p 



E,{\p\cose- E,) 



K 



IpKIpI - Ei,cosi 



K 



(2.36) 



Finally, the transition amplitude can be cast in the form: 



gl\MAJ)\' + Cg, + g,,)M%J)-g.,M.iJ) 
+Coo [3? [{g.MAJ) - 2g,MUJ)) g.Ml.iJ)] + {t, - 2g,g,,) |A<°(J)| 



A/=0,±1 



ig,-Mg^)M^\j)-gyM':^,{J) 



Ml 



+2£0o3? g 



g,Mv{J) - gyMl.{J) 



MUJ) 



+ [gA - gj- 



(2.37) 



{gA + 9pi)M'i{J) - gj^M^iJ) 

The muon capture transition amplitude T\{Jf) can be derived from the result (j2.32|) for 
the neutrino-nucleus reaction amplitude, by keeping in mind that: i) the roles of p and q 
are interchanged within the matrix element of the leptonic current, which brings in a minus 
sign in the last term of C±i^±i, ii) the momentum transfer turns out to he k = q — p, and 
therefore the signs of the right-hand sides of ()2.36p have to be changed, and iii) the threshold 
values {p : k ^ q, k(i, E^ — me) must be used for the lepton traces. All this yields: 

= ^00 = ^90 = 1) -^ii.iti = 1 =F 1- (2.38) 

One should also remember that instead of summing over the initial lepton spins se, as done 
in ()2.20|) . one has now to average over the same quantum number, getting 

A(^/) = §MUJf), (2.39) 

where (pis is the muonic bound state wave function evaluated at the origin, and E^, = 
m^ — (M„ — Mp) — Ej^ — Ef + Ei, where E'^ is the binding energy of the muon in the IS* 
orbit. In view of ()2.39|) the transition amplitude reads: 



UJf) = ^^E[l('^/l|O0(^)-Oo(J)||J.)P + 2|(J^||O^i(J)||J,) 
which, based on the parity considerations, can be expressed as 



(2.40) 



Xv(J/) = G^Y.\\i9y + 9v)My{J)-gyMUJ) + {gA + gj.~g,)M'i{J) - g^M^iJ) 



+ 2 



{gA + gJM/{J)-g^Mj{J) 



(2.41) 
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where gp = g^,^ — g^^. In the case of muon capture, it is convenient to rewrite the effective 



couphng constants as 



9^ 



9^ 



2M 



9 A = 9a 



2M' 



9w — [9v + gm)-^^] 9pi — 9p'^^- 



(2.42) 



Lastly, we mention that the S-values for the Gamow- Teller (GT) beta transitions are 
defined and related to the ft- values as j22i] : 

\9A{Jf\W\\Ji)? 



6146 

B{GT) = ^ sec. 



(2.43) 



2Ji + l ft 

Let us now compare our matrix elements with those currently used in the literature. 
First, in the Walecka's notation (see Eqs. (45.13) in Ref. one has 



O0(J) = Mjo 
Om{J) = 



C 



JO; 



1 

V2 



JM 



for M = 
for M = ±1 



where the meaning of A^jo and Cjo is self evident, while 



Vm = g^SJJM-^'V2 ^'""F 



JLM 



L=J±1 



/-i-mag 



9a^JLM — 9vPjLM 

9a^JJM + 9v^JJM FjlaSjlm- 



The matrix elements defined by Donelly p| are related to ours as: 



L=J±1 



S"?(«:r 



YjA/(Kr), 



K J 



2-^72 (-) ^'''FjLMPLJMiKr) 
'' L=J±l 



L=J±1 



•J-1 



L=J±1 

— YjM(Kr,cr • v). 

K ) 



(2.44) 



(2.45) 



(2.46) 
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The relationship between our formahsm and those from Refs. can be obtained from 



the formula 



(l 1 




II' 1 


'] 


, M 


-M J 


^0 M 


-M j 



3(2L + 1) 



V 2 




+ - 



(: 




:)l 


II 1 


A 














(2.47) 



For the matrix elements of Kuramoto et al. we get 



mm? 



J A/=0,±1 



and for those of Luyten et al. 



(2.48) 



V 



Ml 



Ml 



M/ J 
2 



^) E E 

"V/ J M=0,±1 

— )'ei-^°(^)i'- 



(2.49) 



III. PROJECTED QRPA FORMALISM 

We have shown in Ref. that to account for the weak decay observables in a light 
N = Z nucleus in the framework of the RPA, besides including the BCS correlations, it 
is imperative to perform the particle number projection. It should be remembered that in 
heavy nuclei the neutron excess is usually large, which makes the projection procedure less 
important that in light nuclei 23[. In this section we give a more detailed overview of the 
projected BCS (PECS) and projected QRPA (PQRPA) approximations. 

When the the excited states \Jf) in the final (Z ± 1, =p 1) nuclei are described within 
the PQRPA, the transition amplitudes for the multipole charge-exchange operators Yj, etc, 
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listed in Table ITTl read 

(J/,Z + /i,iV-/x||Yj||0+) 



E 



+ 



K-^{pnJ) 



A^(pnJ) 



with the one-body matrix elements given by 



^Y;{pnJf) 



(3.1) 



A^{pnJ) 



^ip^n, for /i = +1 



V2J + 1 j Mnt>p, for /i = — 1 



(3.2) 



where 



1- 

are the PBCS number projection integrals. 
The PBCS gap equations are 



^ ^ z\lr^''^- a^' = ul + zlvl (3.3) 



2ekUkVk- /^k{u^- Vk) = 



where 



1 _ {2jk' + ly/' 



2 5 (2j, + 1) 



YjjUk'Vk'G{kkk k ; 0) 



are the pairing gaps, and 



Gk = ek 



(3.4) 



(3.5) 



(3.6) 



{2jk + ly/' 

are the dressed single-particle energies. The PBCS correction term Ae^ can be found in Ref. 

and F and G stand for the usual particle-hole (ph) and particle-particle (pp) matrix 
elements, respectively. 

The forward, X^, and backward, Y^, PQRPA amplitudes are obtained by solving the 
RPA equations 



B 



-B* -A! 



L0„ 



Yn 



(3.7) 
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with the PQRPA matrices defined as: 

+ [upUnUp'Un'I^-^^^{pp')l''~^^^{rin') + VpVnVp,Vn'I^-''^^{pp')l'''''^''{rin')]G{pn,pn- J)] , 
B{pn,pn'; J) = N^{pnY^'^N^^{pnY^''^I^~'^{pp)I^~^{nn) 

X [{VpUnUp'Vn' + UpVnVp'Un')F{pn, pri] J) + {UpUnVp'Vn' + VpVnUp'Un')G{pn, pri] J)], (3.8) 

where 

N^ipn) = I^-'+^ipY^'-'+^in), (3.9) 

are the norms, 



(3.10) 



. ^ R^{k) + R^,{kk) _ Bl 
are the projected quasiparticle energies, and the quantities are defined as 

fcl fclfc2 

X [t;^^t;^2F(/ciA;i/c2A;2; 0)I^-\kik2k) + Uk,Vk,Uk2Vk2G{kikik2k2; 0)I^-^{kik2k) 

= Cklull'^ikk) - vll^-\kk)] + E -11^^12 

X {t;2^F(A;iA;iA;A;; 0)[m2/^-2(A;iA;A;) - vll^~\kikk)] - Uk,Vk,UkVkG{kikikk; 0)I^-\kikk)} . 

(3.11) 

It is worth to note that the PQRPA formahsm is vahd not only for the particle-hole charge- 
exchange excitations (Z ± 1, =1= 1), but also for the charge-exchange pairing- vibrations 
(Z ± 1, A^ ± 1). In the later case one has simply to do the replacement fi —>■ —fi in the 
neutron sector. 

The usual gap equations are obtained from Eqs. (j3.4p - (j3.6j) by: 

1. Making the replacement — — Ak, with Ak being the chemical potential, and 
taking the limit ^ 1. That is, the Eq. ()3.4|1 remains as it is, but instead of ()3.5|1 
and (j3.6p one has now 

Afc = E ^^|^^|y^w,,t;fc,G(fcA;fc'A;'; 0), (3.12) 
14 



and 



= e, - Ak + E (2/^ + 1)1/2 vlF{kkk'k'; 0). (3.13) 
2. Imposing the subsidiary conditions 

jp jn 

as the number of particles is not any more a good quantum number. 

Finally, the plain QRPA equations are recovered from ()3.7p and ()3.8|) by: i) dropping the 
index /i and taking the limit — 1, and ii) substituting the unperturbed PBCS energies 
by the BCS energies relative to the Fermi level, i.e. by 

Ef) = ±Efc + Ak, (3.15) 

where Ek = (e| + A|)^/^ are the usual BCS quasiparticle energies. In this way the the 



unperturbed energies in ()3.8|) are replaced as j41| 

4-'^" + ef-'-" - E,^ + E,^ + /i(Ap - A„). (3.16) 

IV. NUMERICAL RESULTS AND DISCUSSION 

In this section, our theoretical results within the PQRPA are confronted with the ex- 
perimental data for the /i^ B C + muon capture rates, as well as for the neu- 
trino cross sections involving the DAR reaction Ve C —>■ ^"^N + e~ and DIF reaction 
i^ti C ^'^N + We also exhibit our predictions for the C N + e~ differ- 

ential cross section for i/e-energies in the DIF energy range. At variance with our previous 
work, we consider here also the velocity dependent matrix elements J^A'iJ) and A^*f(J), 
defined in 



li. 



The calculations shown here were done in the same way as in the previous work 
That is, for the residual interaction we adopted the delta force, 

V = -A7r{v,Ps + VtPt)6{r), (4.1) 



which has been used extensively in the literature |2J, |2^, |2a] for describing the single and 
double beta decays. The configuration space includes the single-particle orbitals with nl = 
(Is, Ip, Id, 2s, If, 2p) for both protons and neutrons. 
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TABLE III: BCS and PBCS results for neutrons. stand for the experimental energies used 

in the fitting procedure, and Cj are the resulting s.p.e. The underlined quasiparticle energies cor- 
respond to single-hole excitations (for jh = lsi/2,lp3/2) a-nd to single-particle excitations (for jp = 
lpi/2, l'^5/2: 2si/2, 1^3/2, 1/7/25 2p3/2, 2^1/2.1/5/2)- ^^'^ nou-underliued energies are mostly two hole-one par- 
ticle and two particle-one hole excitations. The fitted values of the pairing strengths tif"^"" in units of 
MeV-fm^ are also displayed. 



Shell 






BCS 




^3 


PBCS 




1 ,St /o 

^'^1/2 




11.34 


-35.13 


-23.58 


19.93 


-34.99 


-22.37 


lP3/2 


-18.72 


-5.07 


-18.72 


-7.80 


-1.28 


-18.73 


-7.24 


lPl/2 


-4.94 


-4.94 


-18.85 


-2.07 


-4.95 


-22.33 


-1.51 


11^5/2 


-L09 


-1.09 


-22.70 


2.12 


-1.09 


-26.82 


2.16 


2Si/2 


-L85 


-1.86 


-21.93 


2.70 


-1.85 


-25.98 


2.68 


14/2 


2.72 


2.72 


-26.51 


6.24 


2.73 


-30.79 


6.26 


1/7/2 


5.81 


5.82 


-29.61 


8.14 


5.83 


-33.61 


8.17 


2P3/2 


7.17 


7.18 


-30.98 


11.49 


7.16 


-35.23 


11.47 


2pi/2 




12.89 


-36.69 


17.30 


12.89 


-41.01 


17.32 


1/5/2 




16.72 


-40.52 


19.18 


16.72 


-44.58 


19.21 


pair 




23.16 


23.92 



Most of the bare single-particle energies (s.p.e.) ej, as well as the value of the singlet 
strength within the pairing channel {v^°'^^), were fixed from the experimental energies Ej^^ 
of the odd-mass nuclei ^^C, ^^B, ^^C and ^^A^. That is, in the BCS case: 1) we assume that 
the ground states in ^^C and ^^B are pure quasi- hole excitations Ej^^ \ with jh = lp3/2, and 
that the lowest observed 1/2", 5/2+, 1/2+, 3/2+, 7/2~ and 3/2" states in ^^C and ^^A^ are 
pure quasi-particle excitations with jp = lpi/2, lc?5/2, 2si/2, 14/2, 1/7/2) 2^3/2, and 2) 
the s.p.e. of these 7 states and the pairing strength, that appear in the BCS gap equation 
()3.4|) . ()3.12|) - (I3.14|) . were varied in a search in order to account for the experimental 
spectra E^^^ [27[. In the PBCS case we proceed in the same way, i.e. we solve the Eqs. 
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TABLE IV: Same as Table UTTI but for protons. 



Shell 


pexp 

i 


3 


BCS 


^3 


3 


PBCS 

3 


^3 


1 / 2 




14.94 


-33.13 


-21.62 


19.44 


-32.99 


-20.41 


1 Ho In 


-15.95 


-2.24 


-15.96 


-4.95 


1.46 


-15.95 


-4.40 


1 T)^ In 


-1.94 


-1.95 


-16.25 


1.01 


-1.95 


-19.76 


1.56 


Idt; /n 

^^5/2 


1.61 


1.61 


-19.80 


4.76 


1.61 


-23.90 


4.83 


^■31/2 


0.42 


0.42 


—18.61 


4.88 


0.42 


—22.59 


4.88 


^"0/2 


4.95 


4.95 


-23.14 


8.40 


4.95 


-27.37 


8.43 


1/7/2 


8.42 


8.42 


-26.61 


10.72 


8.42 


-30.58 


10.75 


2^3/2 


9.76 


9.77 


-27.96 


14.05 


9.77 


-32.22 


14.07 


2pi/2 




16.23 


-34.43 


20.63 


16.22 


-38.73 


20.66 


1/5/2 




19.70 


-37.90 


22.15 


19.71 


-41.96 


22.20 


pair 
s 




23.13 


23.92 



dH-dSSI), and, instead of fitting and ^ to Ef^ we fit now ej^ and -ef^"^'^"^ We 
have considered the faraway orbitals lsi/2, 2pi/2 and I/5/2 as well, assuming the first one to 
be a pure hole state and the other two pure particle states. Their s^.e. were taken to by 
that of a harmonic oscillator (HO) with standard parameterization 28[. The single-particle 
wave functions were also approximated with that of the HO with the length parameter 
h = 1.67 fm, which corresponds to the estimate hu = 45A^^/^ — 25A~^/^ MeV for the 
oscillator energy. The final results for neutrons are shown in Table UTTI and those for protons 



in Table IIVI It is worth noting that the PBCS neutron-energy ^ip^^^ — -1.28 MeV nicely 
agrees with the experimental energy E^^^- = -1.26 MeV in ^^C. In the same way the PBCS 
)roton-energy ef = 1.46 MeV agrees with the measured energy En in- = 1-55 MeV in ^^N 



. This does not happen in the BCS case, where E[p^^^ = —5.07 MeV for neutrons and 
^ = -2.44 MeV for prot ons, which clearly shows the necessity for the number projection 
procedure. At this point it could be useful to remember that, while in ^^C and ^^B the state 
3/2i is dominantly a hole state, in ^^C and ^^A^ it is basically a two-particle-one- hole state. 
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N = 6 



N = 8 



-30 -20 



-10 O 10 20 30 

e,(MeV) 



FIG. 1: Neutron quasiparticle excitation energies for iV = 6 and N — 8. The states are ordered as lsi/2; 
lP3/2, lPi/2, 2si/2, 1^3/2, 1/7/2, 2p3/2, 2pi/2, and 1/7/2, and the energies are indicated by circles. 



Before proceeding let us remember an important issue in the description of the N = Z 
nuclei within the QRPA, which is more inherent to the model itself that to the occasional 
parameterization that might be employed. In fact, a few years ago Volpe et al. [9] have 
called attention to the inconveniences of applying QRPA to ^^A^, since the lowest state 
turned out not to be the most collective one. Later on we have shown 1J| that the origin of 
this difficulty was the degeneracy among the four lowest proton-neutron two-quasiparticle 
states |lpi/2lp3/2), |lp3/2lp3/2), |lpi/2lpi/2) and 1 1^3/21^1/2)- It also has been shown in Ref. 

that it is imperative to use the projected QRPA for a physically sound description of 
the weak processes among the ground states of the triad \}^B, ^^C, ^^A^}. In fact, when the 
Fermi level is fixed at N = Z = 6, their BCS energies 



J 71 



3p 



Ej^ + Ej^ + Ap - An; for ^^A^ 



3pJn 



El 



(+) 



E. 



Ej„ - Ap + An; for ^'B 



E,^ + E,„ + Ap + An; for '^N 



(4.2) 



-E, 



(-) 



(-) 



Jp Jn 



Xr 



\„: for 



are almost degenerate for all four odd-odd (Z ± 1, A^ =1= 1) and (Z ± 1, A^ ± 1) nuclei. As 
illustrated in Fig. ^ this, in turn, comes from the fact that for N = Z = 6 the quasiparticle 
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energies Eip^^^ and -E'lpg/j ^^^^ very close to each other. In the upper panel of Fig. |21 are 
shown the BCS energies fl4.2|) . as a function of A^, of the just mentioned four states in the 
case of nitrogen isotopes. One notices that for N ^ Z the degeneracy is removed yet only 
partially. But, as seen from the lower panel in Fig. |21 the degeneracy discussed above is 
totally removed when the number projection is done. 




FIG. 2: Unperturbed two-quasiparticle energies £jpj„ for nitrogen isotopes, as a function of N, of the states 
IIP1/2IP3/2), IIP3/2IP3/2), |lPi/2lPi/2) and |lp3/2lpi/2)- In the upper (lower) panel are shown the BCS 
(PBCS) results. 

Moreover, within the PBCS in the case of ^^A^, for instance, we get from numerical 
calculations 

c" — P — F -I- A 

*^lP3/2lP3/2 ~ *^lPl/2lPl/2 ~ *^lPl/2lP3/2 ^' 

^lP3/2lPl/2 - ^lPl/2lP3/2 + (4.3) 

with A = 3.4 MeV. The meaning of this results can be easily disentangled by referring 
to the shell-model and analyzing the particle-hole (ph) limits of the proton-neutron two- 
quasiparticle states, which are pictorially shown in Fig. 01 One sees that, while |lpi/2lp3/2) 
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16.9 
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e = 35.0 (38.5) 34.7 (34.7)^--^ 35.1 (42.1) 34.8 (38.3) 



X 



X X 















34.7 



(34.7) 















-7.3 


( 0.2) 
















E = -7.0 (-3.2) -7.3 ( 0.2) -6.9 (-6.9) -7.2 (-3.5) 



e =11.2(18.7) 10.9(14.7) 11.3(14.8) 11.0(11.0) 



FIG. 3: Schematic representation of the particle-hole limits for the seniority-two pn-states. The zero- 
angular-momentum couplings of two-particles or two- holes are indicateii by a horizontal bracket, ef repre- 
sents the Fermi energy. The unperturbed energies are given in MeV, being the PBCS results given in 
brackets. 

corresponds to a Iplh state in ^^A^, IIP3/2IP3/2) and |lpi/2lpi/2) correspond to 2p2h states, 
and |lp3/2lpi/2) to a 3p3h state in the same nucleus. Therefore one can expect that the 
energy ordering of these states would be given by (j4.3|) . with A = Ajs — Apair being the 
energy difference between the spin-orbit splitting, Ais, and the pairing energy, Apair- A 
similar discussion is pertinent to the remaining three nuclei ^^B, ^'^N, and ^^B. That is, 
one expects that their lowest states would be, respectively, IIP3/2IP3/2), IIP1/2IP1/2) and 
IIP3/2IP1/2), as it happens in the PBCS case but not within the BCS. Finally, it could 
be worthwhile to indicate that the same pn quasiparticle excitation has quite different ph 
compositions in different nuclei. This can be observed by scrutinizing the four columns in 
Fig. El 

The improvement introduced by the PBCS can also be visualized by making a direct 
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E (MeV) 



FIG. 4: The BCS (upper panel) and the PBCS (lower panel) Gamow- Teller strength function in ^'^N . 

calculation of the unperturbed GT strength in ^"^N , given by 

SoAE) ^ 1 E \9a{pM \n) ? o, (4.4) 

The BCS and PBCS results for ^^A^, when folded with t] = 1 MeV, are compared in Fig. 
13] The PBCS energy ordering, given by ()4.3p . is accompanied by the partial shifting of the 
GT strength to higher energies. Therefore it can be said that within the PBCS the GT 
resonance is quenched even at the level of the mean field. 

In view of the above mentioned disadvantages of the standard BCS approach, from now 
on we will mainly discuss the number projection results. In our previous work [lj| we have 
also pointed out that the values of the coupling strengths Vg and Vt within the particle- 
particle (pp) and particle-hole [ph) channels which are used in > Z nuclei (ff^ = f^'**'', 
and ff^ ^ "^s^)) might not be suitable for N = Z nuclei. In fact, the best agreement with 
data in ^^C was obtained when the pp channel is totally switched off, i.e. v^^ = v^^ = 0, 
and three different set of values for the ph coupling strengths jl^, namely: 

Parameterization I {my. vf = v^"" = 24 MeV-fm^ and vf = vf/O.Q = 39.86 MeV-fm^. 
This means that, the singlet ph strength is the same as yP"*'^ obtained from the proton and 
neutron gap equations, while the triplet ph depth is estimated from the relation used by 
Goswami and Pal {2^ in the RPA calculation of ^^C 
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FIG. 5: The results of the PQRPA calculations, as functions of the pp parameter t, are confronted with 



the experimental data taken from Refs. 



1^, 3 3 for: 1 



L) the energy difference uji+ between the 1^ ground 
state in ^^iV and the 0+ ground state in ^^C (upper panel), 2) the S-value for the GT beta transition 
between these two states (middle panel), and 3) the corresponding muon capture rate A(l^) (lower panel). 
The values of v^^ and vf'^ correspond to PII. 

Parameterization II (PII): ff^ = 27 MeV-fm^, and v^^ = 64 MeV-fm^. These values were 
first used in Refs. l2d, no|| and later on in the QRPA calculations of ^^Ca |3,y,|3l|. 

Parameterization III (PHI): v^^ = vf^ = 45 MeV-fm'^. With these coupling constants it 
is possible to reproduce fairly well the energies of the J'^ = Of and li" states in ^^B and 

The results displayed in Fig. El suggest that the choice vf^ = for the pp parameter in 
the 5 = 0, T = 1 channel could be appropriate for the description of the N = Z nuclei. 
They are shown as functions of the parameter 

t * 



pair / \ I pair i \ ' 

Vs \p) + Vs [nj 

together with the experimental data for: 1) the energy difference a;i+ between the l"*" ground 
state in ^^A^ and the O"*" ground state in ^^C, 2) the -B-value for the GT beta transition 
between these two states, and 3) the corresponding exclusive muon capture rate A(lj'") = 
^exc_ rjj^g values of v^^ and wf^ are that from the PII, but quite similar results are obtained 
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with other two sets of parameters. Note that the ph interaction first shifts the energy uji+ 



is increased, we have an opposite attractive effect. That is, the pp interaction diminishes 
uJi+ all up to t ^ 0.6, where the well known collapse of the QRPA approximation occurs. 
At variance with what happens in the case of heavy nuclei, here the values of B(GT) and 
A(15'') basically rise with t, and the agreement with the data is achieved only when the pp 
interaction is totally switched off. Quite generally the nuclear moments ()2.34|) also depend 
weakly on the S = 1, T = channel parameter v^f, for which we adopt as well the null 
value, just to be consistent with our election of ff^. 

Results for the muon capture rates, and the neutrino (z/e,e~) DAR, and (z/^,/i^) DIF 
reaction flux-averaged cross sections are shown, respectively, in Tables |3 Tables IVll and 
Tables En 

The flux-averaged cross section is defined as 



where ^i{Ey) is the normalized neutrino flux. For electron neutrinos this flux was approx- 
imated by the Michel spectrum, and for the muon neutrinos we used that from Ref. j^. 
The energy integration is carried out in the DAR interval nie + ojj^. < A < 52.8 MeV for 
electrons and in the DIF interval + uoj^ < — '^^^ MeV for muons. 

The full PQRPA calculations, which include the relativistic corrections are listed for all 
three parameterizations I, II, and III. On the contrary, the theoretical results involving just 
the velocity-independent operators Yj(/tr) and Sjii^r) are displayed only for the case PII. 
Contributions of the other two operators, PjL^hzr) and Yj(Kr, cr ■ v), to the muon capture 
rates are small (of the order of 5%) as displayed in Table IVl The only exception are the 0~ 
states where the relativistic operator Yo(Kr, cr • v) dominates over the non-relativistic one 
Soi(Kr). We also see from Tables IVTl and IVTTl that the nonlocality effects on the neutrino- 
nucleus reactions are of the order of 1% and therefore they can be neglected. 

The numerical results are sorted according to the order of forbiddeness of the transition 
moments, which can be: allowed (A): J'^ = 0+, 1+, first forbidden (Fl): J'^ = 0~,1~,2~, 
second forbidden (F2): = 2+,3+, third forbidden (F3): J'" = 3~,4~, and so on. The 
response of the three weak processes to successive multipoles is strongly correlated with the 
average momentum transfer, k, involved in each: (z/g, e~)-reaction, K ~ 0.2 fm~^; /^-capture, 



upwards by ~ 1.5 MeV, from its unperturbed value i^ip3/2iP3/2 



16.8 MeV. Then, when t 




(4.5) 
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TABLE V: Experimental and calculated muon capture rate in units of 10"^ s^^. The full PQRPA calcula- 
tions, which include the relativistic corrections are listed for all three parameterizations . Theoretical results 
that involve only the velocity-independent matrix elements are displayed in parentheses in the third column 
for the case PII. The rates are grouped by their degrees of forbiddeness. We show: i) the exclusive rates 
A{J'^), for JJ = 1^ , 1^ ,2^,2^, ii) the multipole decomposition of the rates '^f-A{JJ) for each final state 
with spin and parity JJ , and iii) the inclusive decay rate A'"'^ = X^j" ^(•^7)- column four are listed the 
results of recent shell-model calculations, which are explained in the text. The measured capture rates are 
given in the last column. 



/i capture 




PQRPA 






Shell Model 




Experiment 


rate 


(I) 


(11) 


(III) 


SMU16] 


SM2fl6] 


SM3fl7] 




allowed 




49.3 % 




39.3 % 








A(l^) 


7.52 


6.27(6.50) 


6.27 


11.56 


6.3 


6.0 


6.2 ±0.3 [33j 


E/A(o;) 


3.68 


2.86 (3.15) 


3.77 


0.21 










20.28 18.14 (18.63) 


18.22 


15.43 








1st forbidden 




46.6 % 




55.7 % 








A(ir) 


1.06 


0.49 (0.51) 


0.98 






1.86 


0.62 ±0.20 [34> 35] 


A(2r) 


0.31 


0.18 (0.18) 


0.16 






0.22 


0.18 ±0.10 [M,^ 


E/A(07) 


2.62 


2.35 (0.72) 


2.35 


2.12 








E/A(17) 


11.84 


10.37 (9.51) 


11.37 


12.25 








E/A(27) 


7.78 


7.12 (6.90) 


7.15 


7.79 








2nd forbidden 




3.9 % 




4.6 % 








A(2+) 


0.19 


0.14 (0.16) 


0.15 






0.25 


0.21 ±0.10 [34> 35] 


E/A(2;) 


1.26 


1.09 (0.89) 


1.17 


1.36 








E/A(3+) 


0.63 


0.57 (0.57) 


0.58 


0.46 










48.16 42.56 (40.7) 


44.67 


39.82 


41.9 


33.5 


38 ± 1 



K ~ 0.5 fm~^; and (z/^, /i~) -react ion, k ~ 1 fm~^. As a consequence, in the first case the 
A-moments are by far the dominant ones, contributing with ~ 83.0%, with the remaining 
part of the reaction strength carried almost entirely (~ 16.6%) by the Fl-moments. For 
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TABLE VI: Experimental and calculated flux-averaged cross section for the ^^C{i'e, e^Y'^N DAR reaction 
in units of 10^^^ cm^. The full PQRPA calculations, which include the relativistic corrections, are listed for 
all three parameterizations. Theoretical results that involve only the velocity-independent matrix elements 
are displayed in parentheses for the case PII in the third column. The multipole decomposition E/^e(>/J) 
for each final state with spin and parity , as well as the exclusive, a^^^ = WeiJJ — ^i), and inclusive 
CTg"'^ — J2 J' '^ei'J'f ) cross sections, are shown. In column four are listed the results of a recent shell- model 
calculations, which are explained in the text. The measured cross section are given in the last column. 



(fe, e ) cross 




PQRPA 






Shell Model 




Experiment 


section 


(I) 


(11) 


(III) 


SMU16] 


SM2fl6] 


SM3fl7] 




allowed 




83.0 % 




83.6 % 










9.94 


8.07 (8.00) 


8.17 


20.86 


7.9 


9.3 


8.9 ±0.3 ±0.9 [4] 


E/^e(0+) 


1.92 


1.35 (1.31) 


2.01 


0.00 








E/^e(i;) 


15.98 14.08 (14.22) 


12.14 


22.52 








1st forbidden 




16.6 % 




16.0 % 








E/^e(07) 


0.07 


0.07 (0.05) 


0.07 


0.04 








E/^e(17) 


1.94 


1.59 (1.43) 


1.80 


1.90 








E/^e(27) 


1.66 


1.43 (1.41) 


1.44 


2.36 








2nd forbidden 




0.4 % 




0.4 % 








E/^e(2;) 


0.07 


0.05 (0.04) 


0.05 


0.08 








E/^e(3+) 


0.03 


0.03 (0.03) 


0.03 


0.03 










21.67 18.60 (18.49) 17.54 


26.93 


12.5 


15.1 


13.2 ±0.4 ±0.6 [4] 



/i-capture, the A- and Fl-matrix elements contribute, respectively, with 49.3% and 46.6%, 
while the F2-moments carry only 3.9% of the total transition rate, and the contribution of 
the F3-ones are negligibly small. Finally, the inclusive cross section in the (i/^, /i~)-reaction 
is spread out rather uniformly over several multipoles and it is necessary to include up to 
fourth forbidden moments, with their intensities distributed as follows: 20.0% (A), 39.5% 
(Fl), 27.1% (F2), 10.2% (F3), and 3.1% (F4). The exclusive contributions, coming from the 
ground state 1^^, are quite different in the three cases due to the just mentioned implication 
of the momentum transfer. That is, of the total transition rates for the (i^e, e~)-reaction, 
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TABLE VII: Idem Table EJ but for the averaged exclusive, a^^" = a^{J] = l,^), and inclusive a"^'' = 
J^j" '^tii-Jf): averaged cross sections for the ^^C(i^^, /i")^^Af DIF reaction in units of 10^"*" cm^. 



{iyfj.,H ) cross 




PQRPA 






Sheh Model 




Experiment 


section 


(I) 


(11) 


(III) 


SMl[16j 


SM2[16] 


SM3fl7] 




allowed 




20.0 % 




17.1 % 










0.74 


0.59 (0.56) 


0.59 


1.16 


0.56 


0.9 


0.56 ±0.08 ±0.10 




0.37 


0.26 (0.26) 


0.39 


0.11 










2.68 


2.33 (2.40) 


2.34 


2.95 








1st forbidden 




39.5 % 




36.4 % 










0.03 


0.04 (0.07) 


0.04 


0.07 










3.34 


2.79 (2.84) 


3.10 


3.55 








E/^m(27) 


2.53 


2.28 (2.13) 


2.29 


2.91 








2nd forbidden 




27.1 % 




22.2 % 








e/^m(2;) 


2.55 


2.28 (2.11) 


2.37 


2.59 








E/^M(3t) 


1.34 


1.23 (1.29) 


1.23 


1.39 








3th forbidden 




10.2 % 




17.7 % 








E/ 0^^(37) 


0.74 


0.69 (0.73) 


0.71 


1.77 








E/^m(47) 


0.68 


0.63 (0.63) 


0.63 


1.41 








4th forbidden 




3.1 % 




6.5 % 








E/^m(4) 


0.22 


0.21 (0.20) 


0.21 


0.66 








e/^m(5;) 


0.20 


0.19 (0.19) 


0.19 


0.51 










14.69 12.94 (12.86) 13.51 


17.92 


13.8 


19.2 


10.6 ± 0.30 ± 1.80 



/i-capture, and (z/^, /i^)-reaction, the exclusive contributions are, respectively, 43%, 15%, 
and 5%. 

Let us say a few words on the comparison of our results with the experimental data: 

• /i capture (Table IVj): All exclusive rates A( J'^) with JJ = Ij*", l|f , 2|f , 2j^, are fairly well 
accounted for by the theory, while the inclusive rate, A'"'^, is overpredicted by ~ 10%. 

• (z/g, e~)-reaction (Table ETJ: Even though the exclusive cross section, a^^^, is well 
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reproduced within the PQRPA, the inclusive one, a™*^, is ~ 40% above the data. A 
plausible explanation for this difference could be the fact that we find a very significant 
amount of the GT strength (32%) and the Fermi strength (7%) within the DAR energy 
interval, 18 MeV ^ MeV, where the electron-neutrino flux, ^(.{E^J, changes 

very abruptly, making the inclusive cross section to be very sensitive to the strength 
distribution of low-lying excited states. 

• (z/^, /i~)-reaction fTable IVTT|) : Here also the exclusive cross section, o^^'^, is in full 
agreement with the data, whereas the inclusive one, a™^, is overpredicted by ~ 20%. 

In the last three tables we also confront our PQRPA results with the shell model calcu- 
lations performed by: a) Hayes and Towner [3| , within the model spaces called by them as 
(iii) and (ivj, and which are labeled here, respectively, as SMI and SM2, and b) Auerbach 



and Brown 



17| which we label as SMS. The multipole breakdown in the contributions to the 
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At first glance our 



cross sections from each multipole is only given for the SMI scheme 
results seem to agree fairly well with the shell model ones, particularly when the amounts 
of allowed and forbidden transition strengths are confronted. However, this is not true, as 
can be seen from a more careful analysis of the multipole structure of the transition rates. 

For instance, the total positive parity capture rates with the 1^ state excluded, i.e. 
A(Jj'), are equal to 5.9 and 3.6 in SMI and SMS, respectively, while we get 16.4 
(in units of 10'^ s~^)- Note also that in shell model calculations almost all GT strength is 
exhausted by the ground state transition, while within the PQRPA only 35% of this strength 
goes into the if state. Another important discrepancy is in the Fermi transitions, for which 
we get E/ A(Oj^)= 2.86, while they obtain only E/A(Oj^) = 0.21 (in units of 10^ s'^). 

One sees that our af^^ is consistent with the SM2 and SMS shell model calculations, while 
our a™'^ falls in between the calculations SMI and SM2 and is ~ 20% larger than the SM3 
result. However, same as in the case of muon capture, we find much more ae(l^)-strength 
in the excited states of ^^A^ nucleus than appears in shell model calculations. Namely, we 
obtain that ^e(<//) = 7.44, while in SMI and SMS this quantity is, respectively, 1.77 

and 0.40 (in units of 10~^^ cm^). More, in all three shell model calculations more than 90% 
of the ae(l^)-strength is concentrated in the ^^A^ ground state meanwhile we find only 57%. 

In the same way as in the shell model calculations, within the PQRPA the total DIF cross 
section is mainly built up from forbidden excitations. Nevertheless, although the PQRPA 
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results agree with the SM2 calculation, they are quite small when compared with those 
provided by the SMI and SMS models, for both o^™ and a™'^. The differences in a™'^ come 
not only from the positive parity contributions but also from the negative parity ones. 

As a corollary of the above discussion we would like to stress that it is not easy to asses 
whether the PQRPA results are better or worse than the shell model ones. We just can say 
that with the use of only a few essentially phenomenological parameters the PQRPA is able 
to account for a large number of weak processes in ^^C. (The experimental energies of the 
3/2]^ state in ^^C and ^^N are also predicted within the PBCS.) That is, the s.p.e., ej, for 
j = lp3/2, lpi/2, lc?5/2, 2si/2, lc?3/2, 1/7/2, 2^3/2, and the pairing strength w^"*'' have been fixed 
from of the experimental energies of the odd-mass nuclei ^^C, ^^B, ^^C and ^'^A^, while the 
particle-hole coupling strengths v^^ and vf^ were taken from the previous QRPA calculations 
of ■^^Ca 24, 31 1 . Only the particle-particle couplings v^^ and vf^ has been treated as free 
parameters, and we have used here v^^ = vf^ = 0. It could be interesting to inquire whether, 
this rather extreme parameterization is also appropriate for the description of other light 
N = Z nuclei, such as ^'^A^ and ^^O, which have been discussed recently within a shell model 
scheme ^ 

Finally, we would like to point out that in the DIF neutrino oscillation search 1] an 
excess signal of 

m^'P = 18.1 ±6.6 ±4.0 

events has been observed in the ^^(7(1/6, e~)^^A^ reaction, which are evaluated theoretically 
by the expression 



where 77.3 MeV < A°'^'^ < 217.3 MeV stands for the experimental energy window. The 
oscillation probability reads 

. X 9. .X 9 /l.27Am2L\ 
P.^^.XEu) = sin2(20) sin^ \ j , 

where 9 is the mixing angle between the neutrino mass eigenstates, Am^ is the difference 
in neutrino eigenstate masses squared, in eV^, and L is the distance in meters traveled by 
the neutrino from the source. One sees therefore that the extraction of the permitted values 
of 9 and Am? from experimental data might depend critically on the theoretical estimate 
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FIG. 6: (Color online) Calculated cross section a^^ as function of the neutrino energy. The clashed region 
indicates the experimental energy window. 

of (Je{Ey). So far, the electron cross section obtained within the continuum random phase 
approximation (CRPA) 38| has been used This (JeiEy) is confronted in Figure El with 



our QRPA and PQRPA results calculated with PlI. As can be noticed the PQRPA yields 
a substantially different ae{Ey), inside the experimental energy window for the neutrino 
energy. The consequences of this difference on the confidence regions for sin^(2^^) and Am^ 
will be discussed in a future paper. 

After finishing this work we have learned that quite recently Nieves et al. js^ were able to 
describe rather well the inclusive muon capture rate in ^^C, and the inclusive ^^C(z/^, ^^Y'^N 
and ^^C(z/e, e~Y'^N cross sections, within the framework of a Local Fermi Gas picture which 
includes the RPA correlations. 
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